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Abstract 

The Fermi excitations in hot and dense quark-gluon plasma are studied 
in the Feynman gauge using the temperature Green function technique. We 
find the four well-separated branches for the case m = and establish the 
additional splitting between them (the four different masses) when m 7^ 0. 
The long wavelength limit of these excitations is found in the general case 
of the massive fermions at finite temperature and densities to give the exact 
one-loop spectrum. Simultaneously the many known results are reproduces 
as its different limits. 
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1 Introduction 



At present the interest has essentially grown to study the quark excitation 
spectra (in particular the quark masses) in hot QCD theory at finite fermion 
densities. In the region of the large temperatures and densities normally there 
are no massless particles (all quarks and gluons, at least, have the dynamical 
masses m ~ gT [1]) and these masses being not small influent qualitatively 
on many properties of the quark-gluon medium. Moreover the quarks being 
the Fermi particles are arranged under or near the Fermi sphere and the ap- 
pearance of the diquark pairs is very probably for a rather wide temperature 
region. On the other hand this coupling is not evident at once since the 
particle and hole excitations have the completely different effective masses 
and other properties that is especially pronounced for large m. There are 
a lot paper (starting from [2,3]) which in many aspects discuss the problem 
of the mass generations but today this scenario is well studied only in what 
concerns the dynamical quark mass which can be calculated perturbatively. 
The effective quark mass is usually calculated nonperturbatively to gener- 
ate the chiral phase transition for hot and dense quark-gluon plasma and 
this question is open for discussion. Usually these investigations are rather 
complicated and frequently are accompanied by numerical calculations [4,5,6] 
where many of details are screening. 

The goal of this paper is to flnd analytically the Fermi excitations in hot 
and dense quark-gluon plasma and investigate their peculiarities. We use the 
standard temperature Green function technique and flx the Feynman gauge 
for explicit calculations. The case of a zero damping are only considered and 
we concentrate our attention to solve the fermion dispersion relation in a 
more complete form. We flnd the four well-separated branches for the case 
m — and establish the additional splitting between them (the four different 
masses) when m ^ 0. The long wavelength limit of these excitations is 
investigated in the general case of the massive fermions at finite temperature 
and densities to give the exact one-loop spectrum. We find the dynamical 
quark masses (here only two values) and determine the four effective quark 
masses which define the above mentioned splitting. This splitting show that 
the particle excitations and holes have different masses and that takes place 
for any T, /x- values if bare m ^ 0. Simultaneously a number of known results 
are reproduces as the different limits of the solution found to demonstrate 
its reliability. 
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2 QCD Lagrangian and quark self-energy 



The QCD Lagrangian in covariant gauges has the form 

- fiNj^j,xP + ^{d,Vp' + C^{d,S^' + gr'^V^)d,C' (1) 

where G"^ = df^V^ - dyV^ + gf^^VjlVJ; is the Yang-Mills field strength; 
is a non-Abelian gauge field; ■?/'(and are the quark fields in the SU(N)- 
fundamental representation (|A" are its generators and /"'"^ are the SU(N)- 
structure constants) and C" (and (7") are the ghost Fermi fields. In Eq.(l) n 
and m are the quark chemical potential and bare quark mass, respectively, 
Nf is the number of quarks flavours and a is the gauge fixing parameter 
(a = 1 for the Feynman gauge). The metric is chosen to be Euclidean and 

7^ = 1. 

In 

Our starting point is the exact Schwinger-Dyson equation for the tem- 
perature quark Green function 

G-\q) = G,' + T.{q) (2) 

where the quark self-energy has the standard representation 

^^^^ ^ ^ / (^^'^'^^^ " ?)7.G(p)r.(p, q\p - q) (3) 

Here we calculate S only in the one-loop approximation using the bare Green 
functions in Eq.(3) and fixing the Feynman gauge. The bare functions in 
accordance with Eq.(l) have the form 

r;(p, gb - g) = 7. , G\p) = -'lf_; + ^ , VIM = % (4) 

where p = {{pA + in), p} is the convenient notation. The summation over the 
spinor indices are performed within Eq.(3) using the 7-matrix algebra and 
the result is found to be 

y( , iV^ - 1 r d^p h^p^ + 2m 
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Now we introduce two new functions and rewrite Eq.(5) as follows 

^{q)^ni,K^{q)+mZ{q) (6) 

where K^^q) = q^ a{q) + iUfj, b{q). Here = {1, 0} is the standard medium 
vector and all functions separately depend on 54 an |q|. The decomposition 
(6) is not completed [3,7] but all other functions are not generated in the 
one-loop approximation. Using decomposition (6) we transform Eq.(2) into 
the form 

{q, + K,)^ + m^l + Zy 

which is more convenient futher. Setting the determinant of Eq.(7) to zero, 
we find the dispersion equation 

{q^ + K,y +m2(l + Z)2 = (8) 

which determines the Fermi excitation spectra after the standard analytic 
continuation. 



3 Calculations of the quark self-energy 

Now the summation over the Fermi frequencies p4 = 27rT(n + 1/2) in Eq.(5) 
is performed with the aid of the usual prescription [8] and all terms are 
collected in the convenient form using the simple algebraic trasformations. 
The result is given by 



9\N^ - 1) 
N 



I 



rip bi€p + (i7p + 2m)] 



2(27r)3 \ Lep [q, + + e^) ]^ + - pY 



+ 



n 



P (IpI + - ^54)74 - [?7(q - P) + 2m] 



IPI 
where e 



[q4 + iil^+\p\)? + e, 



p-q 



h.c.{m, jji) ^ —{m, ji) > (9) 



1} ^ and 



n 



± 



v^p^ + m^ is the bare quark energy; = {exp/3|p 
{exp/3( ep ± /x) + 1}^^ are the Bose and Fermi occupation number. 



respectively. 
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Before the last integration over the angles being performed it is desirable 
to separate the functions K{q) and Z{q) to simplify all further calculations. 
Using that TrE(g)/4 = m Z(q) we find the function Z{q) 



N 



(27r)= 



Cp [^4 + ^(/^ + Cp) ]2 + (q - p) 



1 



IpI b4 + ^(/^ + |p| 



+ 



h.c.{ii — > —/i) 



(10) 



p-q 



and then analogously we find the vector Kf^{q). At first, Tr74E(g)/4 = iKi^{q) 
reproduces the ii'4-function 



+1 



nz. 



N 

IpI +/^ 



/ 



2(27r)= 



n 



+ 



[g4 + i(/^ + ep) ]^ + (q-p)^ 



iqi 



h.c.{ii — >• —11) 



IpI [?4 + + IpI) ]^ + eq-p 

and then Tr7„E(g)/4 = iKn{q) gives the vector Kn (n=l,2,3) 



(11) 



Kn{q) 



N 
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dP'p 
2(27r) 



Pn 



ep [q4 + i{l^ + ep) 1^ + (q-p)^ 



nt 



(q-p) 



IpI te4 + i(/^+ IpI) ]2 + eq-p 



+ 



h.c.{ii — > — /i.) 



(12) 



All the further calculations are simple but rather complicated when the 
integration over angles within Eqs.(10)-(12) is performed. Here we use the 
two standard integrals and introduce the new notations to simplify the final 
result. These notations are as follows 



q^ — — {iq^ — /i)^ ± 2ep(ig4 — /i) — 2|p||q| 



q2 _ ^2 _ (^jq^ _ _|_ 2ep(ig4 — /x) + 2|p| |q| 
for integrals with the fermion distribution function and analogously 



(13) 



q^ + — (ig4 


-/X)2±2|p|(^g4 


-/^) -2|p||q| 


q2 _|_ j^2 _ ^^g^ 


-/x)2±2|p|(% 


-/x) + 2|p||q| 



(14) 



and for integrals with boson one. 
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After all algebra is made the function Z{q) has the form 



Z{q4, q) 



9\N^-1) 7d\p\ |p 



N 



I 



47r2 2|q| [ Sp 



— ^ log{aJap) 



nZ — 



n: 



B 



log(— ) - log(aJaB 



(15) 



and the similar form has the function Ki[q) 
iK4{q4, q) = ^-^^^7 j 



N 



T) — 71 



d 


P 




Pl 


87r2 2 


|q| 



log(4) + ^^log(aSaS) 



(16) 



The more complicated calculations are necessary to obtained the vector Kn{q) 
where n = 1,2,3. Here we use a definition Kn{q) = qnK{q) and after that 
the scalar K-function will be calculated. The result has the form 



q^ K{q4, q) = 



+dF ln{apap 



N 



47r2 



1 + 



|p||q| 



|p||q| 



hp In(a^ap) + dp In(^) ^ 



+np IpI [1 + ^1 II . (^B \n{a%aB) + ds ln(- 
^ 8|p||q| ^ a 



B 



|q| 

4|P| 



ln(a^a 



Bj 



(17) 



where hp — 



m 



{iqi — fj,y and dp = 2ep(ig4 — A*) ( the analogously 



hs — <l + m — {iq4 — //) and ds — 2|p|(ig4 — /x)) are the new notations. 



4 The quark excitation spectrum 

Here we use Eq.(8) to find the spectrum of the Fermi excitations of hot 
quark-gluon plasma in the presence of background for the different values of 
the T, ji and m parameters. The cases of a zero damping are only considered 
and due to this fact our analytical continuation is trivial. 
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4.1 The spectrum for all |q| in the case m = 

This is a more simple case which can be exactly considered in the analytic 
form. Here we find the Fermi excitation spectrum for all |g| 7^ within the 
standard high temperature technique keeping /i 7^ 0. To this end we expand 
Eqs.(15)-(17) taking only the leading T^-term with the ///T-corrections and 
solve the dispersion equation selfconsistently exploiting everywhere the new 
mass shell. Our starting point is the expressions for and quantities 
which being equal for the case m = determine all functions which are 
necessary to expand. Using the integration of the quantities over \p\ we 
replace |p| to |p|T and find the following expansions 

ln(a+a-) = + O(^) , In ^ = 2 In |^ + O(i^) (18) 

which are able to simphfy essentially the further calculations. Here ^ = 
{iq4 — /x)/|q| is a convenient variable. 

For our case the dispersion equation (8) can be presented as follows 

[ {iq, - //) - K,]' = (1 + (19) 

and at once to solve it as usual 

(%-//)- ^4 = V\ci\ (l + K) (20) 

where 77 = ±1 and we use the definition K4 = iK^. All functions within 
Eq.(20) should be calculated with the aid of the expansions (18) before we 
solved Eq.(20) explicitly. 

These calculations are easily performed and one finds the following result 

ifta,q) = + + ^B(e-i(l-ai-|^) (21) 

for the function K and analogously for another function 

- K,iq,, q) = In + Is (22) 

Here the integrals are defined to be 

g\N^ - 1) 7d\p\ r n+ + n. 



N 



A/^ - n °° 





g^{N^-l) 7d\p\ n;-n^ 

J 47r2 2 ^^^^ 



N 
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Now one should put the expressions found above into Eq.(20) and perform a 
number of the algebraic transformations to find u) = ^\q\. Here ui = {iq^ — /x) 
and our variable is ^ which is a more convenient then |q|. The latter should 
be excluded with the aid of Eq.(20). The result is the quadratic equation 
with respect to cu 

a;2(e - v[l + Ib B{0] )+u;^Ib + Ik ^'^(0 = (25) 
where the functions A{^) amd B{^) are given by 

^(0=^(l + ^ln|^) 

5(0 = e-^(l-ain|^ (26) 
Keeping the accuracy of calculations our solution of Eq.(25) is found to be 

^(0 = - 4^ ± J/^(^ + ^ln^) (27) 

^ 2(e-r;) 'V ""^e-^ 2 ^ ^ 

which presents the four branches of the Fermi excitations in the medium. 
Here rj = ±1 and we return the physical variable E = ip^. Eq.(27) is a more 
general one-loop result for the case m = and contains all known ones. For 
example, if = 0, the well-known spectrum found in paper [2] 

-iK)=^^-o'(^ + |ln|^) (28) 

is in agreement with Eq.(27) if one puts /i = and squares it. In Eq.(28) 
cj2 = g^T^/Q as it is given by Eq.(23) for iV = 3. Within Eq.(27) 1 < ^ < oo 
and the long wavelength limit corresponds to ^ — > cxo. For this limit one finds 
the very simple formula 




which, indeed, is correct for any T (and for the case T = as well). We 
demonstrate this interesting fact in the next section although, in principle. 
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this situation is known. Eq.(29) represents the more general expression for 
the dynamical quark mass when m = 



N 



47r2 



^ "'v 



(30) 



which for many special cases is well-known. For example, when T = one 
finds the dynamical mass as follows [9] 



TV 167r2 

and in another limit // = this expression was calculated in paper [2]. 



(31) 



4.2 The long wavelength limit with m 7^ 

This is a rather important limit which determines the effective quark mass 
taking into account all radiative corrections: here only the one-loop correc- 
tions. We present our result keeping the perturbative accuracy and establish 
the additional splitting between the branches found above. 

Our starting point is Eq.(8) which for the case m 7^ it is useful to 
rewrite as follows 

iqi{l + a) = /X + 6 + r;/^^^2(7q:^j7:^"^i^7^ (32) 

where 77 = ±1 and all functions should be calculated by using Eq.(6) and its 
decomposition. These calculations yield the simple result 

^(^,^(^^Mi_^l^£H!l^. = - Mate) (33) 

and now Eq.(32) can be transformed in the form 

iq4 = ^ -i{uK{q)) + 77 ^m2(l + Z)^ + + (qK(?)) (34) 

which is convenient for the futher calculations in the case studied. Below 
Eq.(34) will be presented in a more explicit form in the long wavelength 
limit. 
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If the limit |q| = is only considered all functions which determine 
Eq.(34) can be simplify by doing the necessary expansions within Eqs.(15)- 
(17) or exploiting directly Eqs.(10)-(12) in the |q| = hmit. After all al- 
gebraic transformations being made all terms are collected to reproduce the 
final result in the convenient form. 

The function ^4(^4, 0) is found to be 



-i{uK){q^,Q) 



N 



d\v\ 
47r2 



4p^6p 



m 



+ 



2e, 



(1+ ^ 



(35) 



+71, 



B 



4p2 - (% -/X)2[l 



m 



1 + 



(ig4 — — rn? 



2|p| 



(i?4 - IJ)"' 



iq4 - n 



and the function ^(^4, 0) which determines the mass renormalization can be 
presented a follows 



Z{q,,0) 



N 



f d\p\ f 

J 47r2 [ 4el 
^ p 



4p2 



{iq4 - - [{iq4 - /x)^ + m 



212 



+ 2{iq4- fi) 



R 



4|p| [{iqi- f^f 



m 



(36) 



The more complicated calculations yield that (qK)(g4, 0) = 

Now our problem is to solve Eq. (34) explicitly. To this end we should put 
Eq.(35)-(36) on the new mass shell [iq^ — ^) = uo and in this form substitute 
these integrals into Eq.(34). However in this case the arisen equation is very 
complicated and can be used only for numerical calculations [4] . On the other 
hand keeping the perturbative accuracy the obtained integrals are possible 
to simplify considering the quantity uj'^ — ~ inside them. In doing so, 
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the result has the form 



iK4^{q4, 0) 



Ia 



ijj 



- h 



Z(g4,0) = -2Iz + 2-^ 

(jO 



(37) 
(38) 



where the new integrals are: 

g\N^-l)Jd\^\ 



N 



I 





oo 



47r2 



B 



+ IPl 



g'{N'-l) fd\p\ + 



/ 



47r2 2e. 



(39) 



and Ib was earlier determined by Eq.(24). Now the dispersion equation is 
the simple quadratic equation 



u;'^-u;[ri m(l - 2Iz) -Ib]-{Ia + 2r]mlB ) 







(40) 



whose solution reproduces the final result. This result has the form 
1 



r/ m(l - 2Iz) - Ib 



± 



[77 m(l - 2Iz) - Ib 



+ {lA + 27jmlB) (41) 



which gives the additional sphtting of the branches found above. This spht- 
ting show that the particle excitations and holes have different masses and 
that takes place for any T, yU-values if the bare mass m ^ 0. In Eq.(41) all 
parameters are free and this solution extends the previously known result [6] 
taking into account all one-loop corrections. The dynamical quark mass is 
defined to be 



I A + 2-qmlB 



(42) 



and this expression is the more general one-loop result for this quantity. The 
quark mass shift has a rather complicated form 



5m — — - 



77m(l + 27z) + 7b 



± 



'[rim{l-2Iz)-lB, 



+ {I A + 2r]mlB) (43) 
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and is different for eacfi brancfies as it mentioned above. Tliis means tliat a 
question about plienomenological quarlc mass remains open and requires tlie 
additional studying. 

In the limit T = all integrals in Eq.(41) are exactly calculated. For 
example within the standard QCD (where = 3) one has 



61:^ m i-K^ * 

/.=i;(.v^+".M„^i±^^) (44) 

and the spectrum can be investigated by using the simple algebraic transfor- 
mations. In particular case when m = the spectrum has a rather simple 
form 




^ = Ml-£.)±^^ (45) 



which presents the one-loop result for the massless quark excitation in the 
cold medium. In the case m = Eq.(41) coincides with Eq.(29) as it was 
promised. The found equality proves Eq.(29) and this remark is very useful 
for the practical calculations. 



5 Conclusion 

To summarize we have solved the fermion dispersion relation in the general 
case of the massive fermions at finite temperature and densities. We find that 
Fermi excitations in medium for m = have the four well-separated branches: 
two of them present the particle excitations and two other correspond to 
the antiparticles ones. The additional splitting of the branches mentioned 
above is found for the case m 7^ which demonstrate that the effective mass 
for all branches are different for any T, //,- values if only bare m 7^ 0. The 
dynamical mass have only two different values in the general case m ^ and 
an unique value in other cases. Here we see the problem with the definition 
of the phenomenological quark mass since there is a question which branch 
it presents. It is also essential that the particles excitations and the hole 
ones have a different behaviour with respect of the bare mass and there is a 
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problem to find any correlations between them . Another problem concerns 
the gauge dependence of all results found here and in other papers (see 
e.g [10]). This is a serious problem and it is not excluded that a number 
of the one-loop results are, indeed, gauge dependent. In the first place it 
concerns a damping (where resummation is necessary [11]) and probably all 
the ncxt-to-lcading order terms where the same scale ~ g'^T is essential. 
However, in our opinion the long wavelength limit of any excitation spectra 
when the damping is not important should be gauge independent. This 
statement is in agreement with [12] but there is a different result [13] where 
even T^-term is gauge dependent. So today the problem remains and in 
this situation the Feynman gauge used here is a more reliable one to have 
correct result. This problem is still pronounced for many nonperturbative 
calculations which usually use the axial gauge (see e.g. [5, 14]). Moreover at 
present the peculiarities of many singular gauges arc exploited to produce 
the physical results [15]. However these calculations being very complicated 
requires the special attention and should be discussed in a separate paper. 
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